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Abstract 

We analyze circumstances under which the microscopic dynamics of par- 
ticles which are driven by a forced, gradient-type flow can be consistently 
interpreted as a Markovian diffusion process. Special attention is paid to dis- 
criminating between forces that are presumed to act selectively upon diffusing 
particles, while leaving the random medium statistically at rest (Smoluchowski 
diffusion processes), and those perturbing the random medium itself and thus 
creating the nontrivial flows. We focus on the deterministic "stirring" scenar- 
ios. 



To analyze random perturbations that are either superimposed upon or are in- 
trinsic to a driving deterministic motion, quite typically a configuration space equa- 
tion 

X = v{x, t) (1) 

is invoked, which is next replaced by a formal infinitesimal representation of an Ito 
diffusion process 

dX{t) = b{X{t),t)dt + V2DdW{t) . (2) 

Here, W(t) stands for the normalised Wiener noise, and D for a diffusion constant. 

*Presented by P. Garbaczewski at the X Symposium on Statistical Physics, Zakopane, Septem- 
ber 1-10, 1997 
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The dynamical meaning of b{x,t), and thus reasons for making a substantial 
difference between the forward drift of the process and the driving velocity field (1), 
relies on a specific diffusion input and its possible phase-space (e.g. Langevin or 
that coming from deterministic dynamical systems) implementation, that entail a 
detailed functional relationship of v{x,t) and b{x,t), and justify such notions like: 
diffusion in an external force field, diffusion under various strains, diffusion along, 
against or across the driving deterministic flow, [|I|. We shall not touch upon an 
important issue of diffusion under shear, 0, when nontrivial vortices may arise, 
by assuming from the very beginning that only the gradient velocity fields and 
deterministic forces are of interest for us in the present paper. 

The pertinent mathematical formalism corroborates both the Brownian motion 
of a single particle in flows of various origin and the diffusive transport of neutrally 
buoyant components in flows of the hydrodynamic type. However, our major issue 
is a probabilistic interpretation of various linear and nonlinear partial differential 
equations of physical relevance, hence with a slightly abstract flavour put against 
the generally favoured practical reasoning. Expressing that in more physical terms, 
we address an old-fashioned problem of "how much nonlinear", "how much time- 
dependent", and generally-" how much arbitrary" can be the driving velocity field 
to yield a consistent stochastic diffusion process (or the Langevin- type dynamics). 
Another issue is to get hints about a possible non-deterministic origin of such fields, 
i- 

Clearly, in random media that are statistically at rest, diffusion of single tracers 
or dispersion of pollutants are well described by the Fickian outcome of the molecular 
agitation, also in the presence of external force fields (then, in terms of Smoluchowski 
diffusions). On the other hand, it is of fundamental importance to understand how 
flows in a random medium (fluid, as example) affect dispersion. Such velocity fields 
are normally postulated as a priori given agents in the formalism and their (molecular 
or else) origin is disregarded. Moreover, usually the force exerted upon tracers is 
viewed independently from the forcing ("stirring") that might possibly perturb the 
random medium itself and create nontrivial (driving) flows. 

Except for suitable continuity and growth restrictions, necessary to guarantee the 
existence of the process X{t) governed by the Ito stochastic differential equation, 
the choice of the driving velocity field v{x, t) and hence of the related drift b{x, t) is 
normally (in typical physical problems) regarded to be arbitrary (except for being 
"not too nonlinear", see however at van Kampen's discussion of that issue in Ref. 




The situation looks deceivingly simple, [0, if we are (for example) interested in a 
diffusion process interpretation of passive tracers dynamics in the a priori given flow 
whose velocity field is a solution of the nonlinear partial differential equation, be it 
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Euler, Navier-Stokes, Burgers or the like. An implicit assumption, that passively 
buoyant tracers in a fluid have a negligible effect on the flow, looks acceptable 
(basically, in case when the concentration of a passive component in a flow is small). 
Then, one is tempted to view directly the fluid velocity field v{x, t) as the forward 
drift b{x, t) of the process, with the contaminant being diffusively dispersed along 
the streamlines. 

Here, apparent problems arise: irrespectively of a specific physical context and 
the phenomenology (like e.g. the Boltzmann equation with its, as yet, not well 
understood Brownian motion approximation) standing behind the involved partial 
differential equations, some stringent mathematical criterions must be met to justify 
the diffusion process scenario, be it merely a crude approximation of reality. 

Namely, in general, the assumed nonlinear evolution rule for v[x, t) must be 
checked against the dynamics that is allowed to govern the space-time dependence 
of the forward drift field b{x,t) of the pertinent process, |Q, which is not at all 
arbitrary. The latter is ruled by standard consistency conditions that are respected 
by any Markovian diffusion process, and additionally by the rules of the forward 
and backward Ito calculus, Q , the mathematical input that is frequently ignored 
in the physical literature. 

Normally, the pragmatically oriented authors do not pay any attention to such 
problems and feel free to use any (deterministically or not) motivated velocity fields 
as forward drifts. In that case, serious troubles follow. 

Indeed, the closely related issue we have analyzed before, [^, where as a by- 
product of the discussion, the forced Burgers dynamics 

dtVB + {vb ■ V)vb = DAvb + Vn (3) 

and the diffusion-convection equation 

dtc + {vb ■ V)c = DAc (4) 

(originally, for the concentration c{x, t) of a passive component in a flow), in case of 
gradient velocity fields, were found to be generic to a Markovian diffusion process 
input and as generically incompatible with the standard continuity equation in the 
compressible regime. In that case, the dynamics of concentration (in general this 
notion does not coincide with the probability density !) results from the stochas- 
tic process whose density p{x, t) evolves according to the standard Fokker-Planck 
equation 

dtp = DAp-V -{bp) , (5) 
the forward drift solves an evolution equation: 

dtb+{b-V)b=-DAb + Vn , (6) 
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and there holds 



b = VB + 2DVlnp 



(7) 



By combining intuitions which underly the self-diffusion description, [§, with 
those appropriate for probabilistic solutions of the so-called Schrodinger boundary- 
data and next-interpolation problem, 0, 0, §] , the above argument can be general- 
ized to arbitrary conservatively forced diffusion processes, quite irrespectively of a 
physical context in which their usage can be justified. 

Namely, let us consider a density p{x, t) , t > of a stochastic diffusion process, 
solving the Fokker-Planck equation (5). For drifts that are gradient fields, the 
potential VL in Eqs. (3) and (6) [whatever its functional form is), must allow for a 
representation formula, reminiscent of the probabilistic Cameron-Martin-Girsanov 
transformation: 

n(f,t) = 2D[a,$ + i(^ + v-b)] , (8) 

where b{x,t) = 2DV$(x, t). The formula (8) is a trivial identity, if we take for 
granted that all drifts are known from the beginning, like in case of typical Smolu- 
chowski diffusions. Nonetheless, we always end up with a concrete space-time de- 
pendent function Q{x, t) which enters the partial differential equation (6). If we take 
Eq. (6) as a starting point with Q a priori given, its solutions may be then sought 
for in turn (to yield the previous a priori given drifts, if the procedure is consistent). 

Also, the functional properties of Q{x, t) are not an innocent feature of the for- 
malism, since for the existence of the Markovian diffusion process with the forward 
drift b{x,t), we must resort to potentials Q{x,t) that are not completely arbitrary 
functions. Technically, [0], the minimal requirement is that the admissible potential 
is a bounded from below continuous function. This restriction will have profound 
consequences for our further discussion of diffusion in a flow, although nothing seri- 
ous happens if Q is bounded and, for example, is the periodic space-time function. 

Remark 1: If we set p = pi + P2, and demand that p\ ^ p solves the Fokker- 
Planck equation with the very same drift t) as p does, then as a necessary 
consequence of the general formalism, p, 0, the concentration c{x,t) = ^^^^ solves 

an associated diffusion-convection equation dtc + {vb ■ V)c = DAc. Here, the flow 
velocity VB{x,t) coincides with the backward drift of the generic diffusion process 
with the density p{x,t). 

We should clearly discriminate between forces whose effect is a "stirring" of the 
random medium and those acting selectively on diffusing particles, with a negligible 
effect on the medium itself. For example, the traditional Smoluchowski diffusion 
processes in conservative force fields are considered in random media that are sta- 
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tistically at rest. Following the standard (phase-space, Langevin) methodology, let 
us set b{x) = ^K{x), where /3 is a (large) friction coefficient and K represents an 
external Newtonian force per unit of mass ( e.g. an acceleration) that is of gradient 
from, K = —VU. Then, the effective potential Q reads: 

D ^ ^ 

S^-^ + ^V.A- (9) 

and the only distinction between the attractive or repulsive cases can be read out 
from the term V-K. For example, the harmonic attraction/repulsion K = =pax, a > 
would give rise to a harmonic repulsion, if interpreted in terms of VQ, in view 
of n = f^x"^ =F The situation would not change under the incompressibility 

condition (cf. also the probabilistic approaches to the Euler, Navier-Stokes and 
Boltzmann equations, 0). 

Notice that by formally changing a sign of fl we would arrive at the attractive 
variant of the problem, which is however incompatible with the diffusion process 
scenario in view of the unboundedness of —Q from below. 

We have thus arrived at the major point of our discussion: a priori, there is no 
way to incorporate the attractive forces which affect (drive) the flow and nonetheless 
generate a consistent diffusion- in- a- flow transport. 

Clerly, there is no reason to exclude the attractive variants of the potential Q 
from considerations, since the deterministic motion is consistent with them. 

Concluding, if the diffusion is to be involved we need to save the situation some- 
how, and this can be made only by incorporating the hitherto not considered "pres- 
sure" term effects. That is suggested by the general form of the compressible Euler 
{F = —W stands for external volume forces and p for the fluid density that itself 
undergoes a stochastic diffusion process): 

dtVE + {ve ■ V)ve = F- -VP (10) 

P 

or the incompressible, [0, Navier-Stokes equation: 

dtVNs + {vNs ■ ^)VNS = -/^vns + F- -VP , (11) 

p p 

both to be compared with the equations (1) and (4), that set dynamical constraints 
for respectively backward and forward drifts of a Markovian diffusion process ? 

Let us stress again that the acceleration term F in equations (10) and (11) 
normally is regarded as arbitrary, while the corresponding term Vfi in (3), (6) and 
(8) involves a bounded from below function. 
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Since, in case of gradient velocity fields, the dissipation term in the incompress- 
ible Navier-Stokes equation (11) identically vanishes, we should concentrate on an- 
alyzing the possible "forward drift of the Markovian process" meaning of the Euler 
flow with the velocity field ve, (10). 

At this point it is useful, at least on the formal grouds, to invoke the standard 
phase-space argument that is valid for a Markovian diffusion process taking place in 
a given flow v{x^ t) with as yet unspecified dynamics nor physical origin. We account 
for an explicit force exerted upon diffusing particles, while not necessarily directly 
affecting the driving flow itself. Namely, P, |], let us set for infinitesimal increments 
of phase space random variables: 

dX{t) = V{t)dt 

dV{t) = p[v{x, t) - V{t)]dt + K{x)dt + /3V2DdW{t) . (12) 

Following the leading idea of the Smoluchowski approximation, we assume that 
(3 is large, and consider the process for times significantly exceeding Then, an 
appropriate choice of the velocity field v{x, t) (boundedness and growth restrictions 
are involved) may in principle guarantee, the convergence of the spatial part 
X{t) of the process (12) to the Ito diffusion process with infinitesimal increments 
(where the force K effects can be safely ignored if we are interested mostly in the 
driving motion): 

dX{t) = v{x, t)dt + V2DdW{t) . (13) 

However, one cannot blindly insert in the place of the forward drift v{x, t) any 
of the previously considered bulk velocity fields, without going into apparent con- 
tradictions. Specifically, the equation (4) with w(x, t) ^ b{x, t) must be valid. 

By resorting to velocity fields v{x, t) which obey Av{x, t) = 0, we may pass from 
(6) to an equation of the Euler form, (10), provided (8) holds true and then the 
right-hand-side of (6) involves a bounded from below effective potential Q. 

An additional requirement is that 

F- -VP = Vn . (14) 
P 

Clearly, in case of a constant pressure we are left with the dynamical constraint 
(6 ^ Ve): 

dtb + ib-V)b = F = Vn (15) 

combining simultaneously the Eulerian fluid and the Markov diffusion process in- 
puts, if and only if F is repulsive, e.g. —V{x,t) is bounded from below. Quite 
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analogously, by setting F = 0, we would get a constraint on the admissible pressure 
term, in view of: 

dtb+(b-V)b=--VP = Vn . (16) 
P 

Both, in cases (15), (16) the effective potential fl must respect the functional 
dependence (on a forward drift and its potential) prescription (8). In addition, the 
Fokker-Planck equation (5) with the forward drift VE{x,t) = b{x,t) must be valid 
for the density p{x,t). 

To our knowledge, in the literature there is only one known specific class of 
Markovian diffusion processes that would render the right-hand-side of Eq. (10) 
repulsive but nevertheless account for the troublesome Newtonian accelerations, e.g. 
those of the from — V^, with +V bounded from below. Such processes have forward 
drifts that for each suitable, bounded from below function V{x) solve the nonlinear 
partial differential equation: 

dtb+{b-V)b = -DAb + V{2Q -V) (17) 

with the compensating pressure term: 

Q = 2^2-^ = ^^ + DV-u (18) 

—t 

u(x, t) = DVln p{x, t) 

Their discussion can be found in Refs. 0, |^, 0, 0. 
Clearly, we have: 

F= -W, V2Q = --VP (19) 

P 

where: 

P(f , t) = -2D^p{x, t) A In p{x, t) (20) 

Effectively, P is here defined up to a time-dependent constant. Another admissible 
form of the pressure term reads (summation convention is implicit): 

^Vk[p{2D^d,dk)lnp] = V,(2Q) (21) 

If we consider a subclass of processes for which the dissipation term identically 
vanishes ( a number of examples can be found in Refs. 0): 

Ab{x, t) = (22) 



7 



the equation (17) takes a conspicuous Euler form (10), ve ^ b. 

Let us notice that (20), (21) provide for a generahsation of the more famihar, 
thermodynamically motivated and suited for ideal gases and fluids, equation of state 
P ~ /9. In case of density fields for which —Aln p ~ const, the standard relationship 
between the pressure and the density is reproduced. In case of density fields obeying 
—Aln p = 0, we are left with at most purely time dependent or a constant pressure. 
Pressure profiles may be highly complex for arbitrarily chosen initial density and/or 
the flow velocity fields. 

To conclude the present discussion let us invoke Refs. H, . The problem of a 
diffusion process interpretation of various partial differential equations is known to 
extend beyond the original parabolic equations setting, to general nonlinear velocity 
field equations. On the other hand, the nonlinear Markov jump processes associated 
with the Boltzmann equation, in the hydrodynamic limit, are believed to imply 
either an ordinary deterministic dynamics with the velocity field solving the Euler 
equation, or a diffusion process whose drift is a solution of the incompressible Navier- 
Stokes equation (in general, without our curlv = restriction), H, ||. The case of 
arbitrary external forcing has never been satisfactorily solved. 

Our reasoning went otherwise. We asked for the admissible space-time depen- 
dence of general velocity fields that are to play the role of forward drifts of Markovian 
diffusion processes, and at the same time can be met in physically signicant con- 
texts. Therefore various forms of the Fokker-Planck equation for tracers driven by 
familiar compressible velocity fields were discussed. 

Our finding is that solutions of the compressible Euler equation are appropriate 
for the description of the general non-deterministic (e.g. random and Markovian) 
dynamics running under the influence of both attractive and repulsive stirring forces, 
and refer to a class of Markovian diffusion processes orginally introduced by E. 
Nelson, |^, Q . That involves only the gradient velocity fields (a couple of issues 
concerning the curl b velocity fields and their nonconservative forcing have been 
raised in Refs. |^). 

Remark 2: Let us stress that a standard justification of the hydrodynamic limit 
for a tracer particle invokes a Brownian particle in an equilibrium fluid. An issue of 
how much the tracer particle disturbs the fluid (random medium) locally and how 
far away from the tracer particle the thermal equlibrium conditions regain their 
validity, normally is disregarded. Moreover, in the standard derivation of local 
conservation laws from the Boltzmann equation, the forcing term on the right-hand- 
side of the Euler or Navier-Stokes equation up to scalings does coincide with the 
force acting on each single particle comprising the system. Thus, in this framework, 
there is no room for any discrimination between forces acting upon tagged particles 
and those perturbing the spatial flows (once on the level of local averages). 
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Quite on the contrary, the force term in the Kramers equation and this appear- 
ing in the related local conservation law for the forward drift or for the current 
velocity of the diffusion process are known not to coincide in general. Typically, the 
action of an external force is confined to diffusing (tagged) particles with no global 
or local effect on the surrounding random medium, cf. standard derivations of the 
Smoluchowski equation. This feature underlies problems with the diffusion process 
interpretation of general partial differential equations governing physically relevant 
velocity fields. Specifically, any external intervention (forcing) upon a stochastically 
evolving (in the diffusion process approximation) system gives rise to a perturbation 
of local flows, which seldom can be analyzed as forcing of any definite type on the 
molecular level. The Smoluchowski theory is here a notable exception, but there 
one has no room for genuine flows and velocity field profiles which are generated in 
the random medium. 



Remark 3: It seems worthwhile to mention a close connection of the considered 
framework with the general issue of executing small random perturbations on the 
level of the classical Hamilton- Jacobi dynamics, , with the related issue of an opti- 
mal control of stochastic processes and with that of the so-called "viscosity solutions" 



of the Hamilton- Jacobi equation, [|T0|. In fact, our (Feynman-Kac, see 0) potentials 
(8), (9) were introduced on the basis of probabilistic arguments via the Girsanov 
or Cameron-Martin theorems about transformations of drifts of the diffusion pro- 
cess. However, an implicit assumption that drifts are defined in terms of gradients 
of suitable logarithmic functions: h = 2DVlog9 and vb = K = —2DVlog9^. (here, 
we employ the notation of our previous publications, [^, where 9^ is a bounded 
solution of the forward generalized diffusion equation, while 9 that of its time ad- 
joint) implies that the compatibility condition (8) can be rewritten in two equivalent 
forms, both involving the modified Hamilton- Jacobi equations. Namely, let us set 
2Dlog9 = $ and -2Dlog9^ = Then, we have: = + ||V$|2 + DA$ 

and at the same time Q = dt<^^ + i|V$*p - DA^^. The latter one is identified 
as the so-called Hamilton- Jacobi-Bellmann programming equation in the optimal 
control of stochastic diffusion processes, [^0|, ^, and via the Hopf-Cole logarithmic 
transformation (take the gradient) is linked to the Burgers equation (3). An issue of 
viscosity solutions of the standard Hamilton- Jacobi equation has been extensively 
studied in the literature as the D I limit of solutions of the modified (e,.g. Bell- 
mann) equation. It is thus clear, on the basis of our previous discussion, that an 
apparent obstacle is hidden in the assumption that a diffusion process is involved. 
Then, suitable restrictions upon Q must be respected, and the attractive versus 
repulsive potential problem appears. 
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